UNCLASSIFIED 

«g|  4  5  0  2  0  5 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION  ALEXANDRIA,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  vith  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


TALOGED  BY  DDC  45020 


|  CONTRACT  TECHNICAL  NOTE 

v 


A  NUMERICAL  INVESTIGATION  OF  THE  VALIDITY 
OF  THE  BORN  APPROXIMATION  FOR 

^TERMINATION  OF  THE  REFLECTION  COEFFICIENTS 
OF  UNDERDENSE  PLASMA  SLABS 


CONTRACT  NO.  DA.04-495-ORD-3567(Z) 
HYPERVELOCITY  RANGE  RESEARCH  PROGRAM 
A  PART  OF  PROJECT  "DEFENDER" 


GM  DEFENSE  RESEARCH  LABORATORIES 


SANTA  BARBARA,  CALIFORNIA 


AEROSPACE  OPERATIONS  DEPARTMENT 


CTN64-08 


,0@/  ZiskK 

SGJpFf 

9Wair^ 


OCTOBER  1964 


NOTICE 


The  information  presented  in  this  Technical  Note  is  pre¬ 
liminary  and  should  not  be  interpreted  as  representing  the 
final  opinion  of  GM  Defense  Research  Laboratories  or  the 
author  on  the  subject.  Further,  to  expedite  prompt  distri¬ 
bution  of  the  technical  contents,  no  attempt  at  an  editorial 
review  has  been  made. 


Copy  No  IPS, 


ARPA  ORDER  NO.  357-63 
PROJECT  CODE  NO.  7400 

S  CONTRACT  TECHNICAL  NOTE  J 


A  NUMERICAL  INVESTIGATION  OF  THE  VALIDITY 
OF  THE  BORN  APPROXIMATION  FOR 
DETERMINATION  OF  THE  REFLECTION  COEFFICIENTS 
OF  UNDERDENSE  PLASMA  SLABS 

S.  ZIVANOVIC 


THIS  RESEARCH  WAS  SUPPORTED  BY  THE 
ADVANCED  RESEARCH  PROJECTS  AGENCY 
DEPARTMENT  OF  DEFENSE 
AND  WAS  MONITORED  BY  THE 
U.S.  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL,  ALABAMA 


GM  DEFENSE  RESEARCH  LABORATORIES 

SANTA  p/ JIBARA,  CALIFORNIA 

aerospace  operations  debartment 


CONTRACT  NO.  DA-04-495-ORD.3567(Z) 
HYPERVELOCITY  RANGE  RESEARCH  PROGRAM 


A  PART  OF  PROJECT  "DEFENDER 


•  • 


CTN64-08 


DDC  AVAILABILITY  NOTICE 

Qualified  requesters  may  obtain 
copies  of  this  report  from  DDC 


OCTOBER  1964 


ABSTRACT 

The  validity  of  the  Born  approximation  for  determination  of 
the  reflection  coefficients  of  plasma  slabs  is  numerically 
examined  for  two  cases:  the  uniform  underdense  plasma  slab 
and  a  plasma  slab  with  a  sinusodial  electron  density  perturb¬ 
ation  superimposed  on  a  constant  electron  density.  The  regions 
of  validity  of  the  Born  approximation  are  discussed  in  each  case. 


i 


CTN64-08 


CONTENTS 

Page 

ABSTRACT  i 

I.  INTRODUCTION  1 

n.  THEORY  2 

Reflection  and  Transmission  of  EM  Waves  from 
Nonuniform  Underdense  Plasma  Slabs  2 

m.  DISCUSSION  OF  RESULTS  4 

IV.  CONCLUSIONS  10 

REFERENCES  10 

ILLUSTRATIONS 

Figure  No.  Title 

1  Geometry  of  the  Plasma  Slab  2 

2  Percentage  Error  in  Reflection  Coefficient  of 
Uniform  Plasma  Slab  for  Born  Approximation 
as  a  Function  of  Normalized  Plasma  Frequency, 

Normalized  Thickness  of  the  Slab,  and  Normalized 
Collision  Frequency  0,  =  0  6 

3  Percentage  Error  in  Reflection  Coefficient  of 
Uniform  Plasma  Slab  for  Born  Approximation 
as  a  Function  of  Normalized  Plasma  Frequency, 

Normalized  Thickness  of  the  Slab,  and  Normalized 
Collision  Frequency  fi,  =  0. 1  7 

4  Percentage  Error  in  Reflection  Coefficient  of 
Uniform  Plasma  Slab  for  Born  Approximation  as 
a  Function  of  Normalized  Plasma  Frequency, 

Normalized  Thickness  of  the  Slab,  and  Normalized 
Collision  Frequency  =  0.  2  8 

5  Reflection  Coefficient  (Born  Approximation)  of  a 
Plasma  Slab  with  Harmonic  Density  Variation  as 
a  Function  of  the  Normalized  Density  Wavelength 

and  a  Density  Modulation  h  =  0.  05  H 

6  Reflection  Coefficient  (Born  Approximation)  of  a 
Plasma  Slab  with  Harmonic  Density  Variation 

as  a  Function  of  the  Normalized  Density  Wavelength 
and  a  Density  Modulation  h  =  0. 1  12 


ii 


CTN64-08 


ILLUSTRATIONS  (Cont'd) 

Figure  No.  Title  Page 

7  Reflection  Coefficient  of  a  Plasma  Slab  with 
Harmonic  Density  Variation  as  a  Function  of 
the  Normalized  Density  Wavelength  and  a 
Density  Modulation  h  =  0.  2 

Dashed  Line  -  Exact  Values  of  Reflection 
Coefficient 

Solid  Line  Born  Approximation  of  the 

Reflection  Coefficient  13 

8  Reflection  Coefficient  (Born  Approximation) 
of  a  Plasma  Slab  with  Harmonic  Density 
Variation  as  a  Function  of  the  Normalized 

Density  Wavelength  and  a  Density  Modulation  h=  0.  4  14 

9  Reflection  Coefficient  (Born  Approximation)  of 

a  Plasma  Slab  with  Harmonic  Density  Variation 
as  a  Function  of  the  Normalized  Density  Wave¬ 
length  and  a  Density  Modulation  h  =  0. 6  15 

10  Reflection  Coefficient  of  a  Plasma  Slab  with 
Harmonic  Density  Variation  as  a  Function  of 
the  Normalized  Density  Wavelength  and  a 
Density  Modulation  h  =  1.  0  and  a  =  5. 25 

Colored  Curve:  Exact  Values  of  Reflection 
Coefficient 

Black  Curve:  Born  Approximation  of  the 

Reflection  Coefficient  16 

11  Reflection  Coefficient  of  a  Plasma  Slab  with 
Harmonic  Density  Variation  as  a  Function  of  the 
Normalized  Density  Wavelength  and  a  Density 
Modulation  h  =  1. 0  and  a  =  5. 5 

Colored  Curve:  Exact  Values  of  Reflection 
Coefficient 

Black  Curve:  Born  Approximation  of  the  17 

Reflection  Coefficient 


iii 


CTN64-08 


I.  INTRODUCTION 

In  the  study  of  the  propagation  of  electromagnetic  waves  through  ionized  gases 
it  is  seldom  possible  to  get  analytic  solutions  which  are  easy  to  interprete 
without  tedious  computing.  Direct  numerical  solutions  give  accurate  results, 
but  they  also  fail  to  give  physical  insight  into  the  mechanism  of  propagation. 
However,  sometimes  it  is  desirable  to  have  an  analytic  solution  which  is  easy 
to  follow  even  if  some  accuracy  has  to  be  sacrificed. 

The  Born  approximation,  applied  to  the  problem  of  EM  scattering  from  an 

underdense  plasma,  serves  this  purpose.  To  use  the  Born  approximation  one 

has  to  assume  that  the  incident  wave  is  not  changed  by  the  medium  through 

which  it  passes,  and  to  neglect  multiple  scattering.  It  has  been  used  often 

(1)* 

and,  for  example,  Booker  and  Gordon'  'have  used  it  successfully  to  solve  the 
problem  of  scattering  from  small  random  inhomogeneities  in  the  troposphere. 

A  good  discussion  of  the  Born  approximation  is  given  in  Schiff^and  other 
textbooks  on  quantum  mechanics  and  wave  propagation. 

In  this  note,  the  Born  approximation  is  used  to  derive  the  transmission  and 
reflection  coefficients  of  uniform  and  nonuniform  plasma  slabs  for  normal 
incidence.  K  is  shown  that  the  Born  approximation  gives  very  good  results 
for  the  reflection  coefficient  of  uniform  slabs  (within  8%  up  to  plasma  frequencies 
of  0.  2  of  incident  frequency  and  slab  thicknesses  up  to  five  wavelengths  of 
incident  wave).  The  transmission  coefficient,  calculated  from  the  Born  approxi¬ 
mation,  has  a  magnitude  greater  than  unity  but  its  phase  shift  is  approximately 
correct.  The  reflection  coefficient  for  a  plasma  in  which  there  are  harmonic 
variations  in  electron  density  is  also  calculated  and  compared  with  the  exact 
results. 


*  Raised  numbers  in  parentheses  indicate  references  listed  at  the  end  of 
this  report. 
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In  this  latter  case  the  errors  in  the  Born  approximation  are  due  mostly  to 
the  scattering  from  the  outside  boundaries  of  the  plasma.  It  should  be 
mentioned  that  a  strong  resonance  occurs  in  the  case  when  the  wavelength  of 
the  fluctuations  is  exactly  equal  to  one-half  of  the  free-space  wavelength. 
This  result  can  be  useful  in  studying  the  reflection  from  the  turbulent  wakes 
of  hypersonic  projectiles  -  a  subject  which  will  be  treated  in  a  subsequent 
note. 


II.  THEORY 

Reflection  and  Transmission  of  EM 
Waves  from  Nonuniform  Underdense 
Plasma  Slabs 

Let  a  plasma  be  confined  to  a  slab  of 
thickness  d,  located  such  that  the  z- 
axis  is  perpendicular  to  the  slab,  as 
depicted  in  Figure  1.  Let  a  plane 
transverse  electromagnetic  wave 
propagating  in  the  direction  of  the 
positive  z-axis  with  its  electric 
field  E.  in  the  x  direction  be 

i 

incident  on  the  plasma  slab.  The 
electric  field  is  given  by 

Ei  =  lx  Eq  exp  [j(wt-kz)] 


(1) 


where  Eq  is  the  complex  amplitude  of  the  incident  wave,  to  is  the  angular 
frequency  of  the  wave,  t  is  the  time,  k=2irX  *  and  X  is  the  wavelength  of  the 
wave  in  free  space. 


It  will  be  assumed  that  the  first  Born  approximation  is  valid,  i.  e. ,  that  the 
variations  of  E.  can  be  neglected  throughout  the  plasma  and  also  that  the 
multiple  scattering  can  be  neglected,  hence  the  incident  field  is  the  total  field 
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throughout  the  plasma.  It  will  also  be  assumed  that  the  properties  of  the 
plasma  change  only  in  the  z -direction. 


A  uniform  layer  of  plasma  of  thickness  a  will  have  the  reflection  coefficient^ 

Er  (1-f/f  )[  1-exp  (-2jk  yJ77T  a)] 

R  ~  E  ~  i -  ~o - - — - 5 - - - -  (2) 

1  (V^o+1>  "  !)  exp(-2jk^77c0  a) 

where  Er  and  E.  are  the  complex  amplitudes  of  the  reflected  and  incident 
electromagnetic  waves,  respectively,  evaluated  at  the  incident  side  of  the 
slab,  and  e  and  are  the  permittivity  of  the  plasma  and  free  space,  respectively. 


For  an  elementary  slab  of  thickness  df  ,  located  at  f,  expression  (2)  reduces 
to  (subscript  "e"  stands  for  "elementary") 


R 

e 


jkdf  =  -^Adf 


where 


(3) 

(4) 


or,  of  one  defines  the  reflection  coefficient  as  the  ratio  of  the  reflected-wave 
amplitude  to  the  transmitted -wave  amplitude  at  the  origin, 


df 


(5) 


Within  the  validity  of  the  Born  approximation  the  reflection  coefficient  of  the 
entire  slab  is 


f  Ae-2*fdr 

A=o 


(6) 
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The  transmission  coefficient  of  a  uniform  plasma  slab  is 
_  4Ve7T  exp  \-  jk  (\fe7e  -  l)al 

T  =  —  ~  — - - - != _  0  J _  /m 

Ei  (yje / + 1)2  -  (-y/cTT  -  l)2  exp  J^-2jk y/e/ eQ  a  j 

where  E^  is  the  amplitudi  of  the  transmitted  wave  and  Et  and  E.  are  evaluated 
at  the  same  point.  For  infinitesimal  thickness 

Tg  =  1-jk  d^[(>/777o-  1)  -  I  (y/777-  l)2]  (8) 

Retaining  only  the  first  power  in  ( yje /e  - 1) ,  (8)  becomes 

Te  =  1_;ik(  V<7T- i)  dr  (9) 

The  change  in  the  transmitted  wave  due  to  the  elementary  plasma  slab  is 
represented  by  the  second  term  in  the  right-hand  side  in  (9).  The  total 
transmission  coefficient  is  hence 


Te  =  Hk  /  (V*/V1)df 


The  phase  shift  in  the  transmitted  wave  is  hence  (for  negligible  losses,  that 
is,  pure  real  e  ) 


Zl 


-  arc  tan  k 


/<V^o- 


1)  df  »  -  arc  tan  £  /  AdT 


*/' 


m.  DISCUSSION  OF  RESULTS 

The  magnitude  of  the  transmission  coefficient  for  real  e  in  the  Born  approxi¬ 
mation  is  always  greater  than  one,  and,  hence,  physically  meaningless. 
However,  for  small  A  (10)  will  give  the  phase  shift  with  an  accuracy  which 
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will  depend  on  the  thickness  of  the  medium  and  the  magnitude  of  A  .  Also 
(10)  gives  an  indication  that,  at  least  for  small  A,  the  angle  of  T  depends 
only  on  the  average  properties  of  the  medium. 


The  reflection  coefficient  (6)  is  a  function  of  the  changes  of  A  throughout  the 
plasma  and,  for  small  A,  the  errors  made  by  using  (6)  instead  of  the  exact 
reflection  coefficient  are  small.  To  illustrate  this  point,  on  Figures  2  to  4, 
the  error 


e  = 


R-R 

a 

R 


100 


(12) 


where  R  is  the  exact  value  of  the  reflection  coefficient  for  the  uniform 
plasma  and  R&  is  computed  from  (6)  assuming  constant  A,  is  plotted.  In  a 
slightly  ionized  plasma^ 

O2 

A='HiT  (13) 

2  2  2 

where  0^  =(q  n)/  (eQ  mw  ),  q  is  the  charge  of  an  electron,  m  is  its  mass, 

n  is  the  electron  density,  and  vc  is  the  collision  frequency 

between  electrons  and  neutral  atoms.  It  is  seen  that  the  error  does  not  exceed 

a  few  percent  for  12  0.  2  and  f2  ^  0.  2  . 

p  c 


If  the  electron  density  has  a  periodic  component,  i.  e. ,  it  can  be  expressed 
as 

n  =  nQ  +  nj  sin  (kjZ+9)  (14) 

where  nQ  >  ^  and  n ^  and  kj  are  constants,  then,  assuming  constant 
collision  frequency,  A  can  be  expressed  as 

n  2 

A=  “i4^“(1+h  sinkiz)  (15) 
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2  2  2 

where  0  "'•-(q  no)(fQmw)  and  h  =  n^/n  .  Substituting  expression 
(15)  into  (6)  gives 


2  r 

R  =  /  (1 +h  sin  )  e  Jkdf 

2i(l-jfic)f-o  1 


(16) 


Integrating,  one  has 


0  ' 
R  =— E2_ 


4(1-  JOc) 


l-e”^2ka  +h(ka) 


-eja(k1-2k)_1  e-ja(k1  +  2k)_1-, 

.ja(k1-2k)  ja(k1 +2k)  J 


(17) 


The  first  term  on  the  right-hand  side  represents  the  Born  approximation  for 
the  reflection  from  a  uniform  plasma,  i.  e. ,  for  h=o.  The  second  term 
represents  the  contribution  due  to  the  periodicity  of  the  plasma. 

It  is  obvious  that  for  kj“2k  a  strong  resonance  occurs.  For  kj*2k*4ffA 
expression  (18)  becomes 


and  hence,  a  ripple  in  electron  density  equal  to  only  20%  of  the  average 
electron  density  can  enhance  the  reflection  of  the  plasma  ten  wavelengths 
long  by  over  six  times  (about  16  db  in  power). 


Off  resonance,  the  contribution  of  the  harmonic  part  of  the  electron  density 
decreases  rapidly  and  fluctuates  around  the  reflection  coefficient  for  a 
uniform  plasma. 
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To  illustrate  the  accuracy  of  the  Born  approximation,  in  Figures  5  to  11 
the  square  of  the  magnitude  of  the  reflection  coefficient  for  the  plasma  with 
electron  density  distribution  given  by  (14)  is  plotted  versus  the  density  wave¬ 
length  M  =  2ir/kj  ,  as  calculated  from  expression  (17).  Then,  for  h  =  0.  2 
and  h  =  1  the  exact  values  of  reflection  coefficient  obtained  by  using  an 
exact  numerical  technique  'are  given  for  comparison.  It  is  seen  that  the 
agreement  between  the  two  curves  is  good  and  that  the  main  difference  comes 
from  the  uniform  component  of  electron  density  (change  of  the  effective  thick¬ 
ness  of  plasma).  Around  the  resonances  (M=  1/2) agreement  is,  for  most 
practical  purposes,  very  good.  The  vertical  scale  is  expressed  in  decibels 
below  the  unity  reflection  coefficient  because  the  square  of  the  reflection 
coefficient  as  defined  by  (2)  is  the  power  reflection  coefficient. 

IV.  CONCLUSIONS 

It  is  seen  that,  depending  on  the  required  accuracy,  the  Born  approximation 
can  be  a  very  satisfactory  approach  for  studying  the  reflections  from  under- 
dense  plasmas  -  up  to  plasma  frequencies  of  about  0.  2  of  incident  frequency. 

Its  chief  merit  lies  in.  its  simplicity  and  ability  to  handle  nonuniform  plasmas. 
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Figure  5  Reflection  Coefficient  (Born  Approximation)  of  a  Plasma  Slab  with 

Harmonic  Density  Variation  as  a  Function  of  the  Normalized  Density 
Wavelength  and  a  Density  Modulation  h  =  0.  05 


Figure  8  Reflection  Coefficient  (Born  Approximation)  of  a  Plasma  Slab  with  Harmonic 
Density  Variation  as  a  Function  of  the  Normalized  Density  Wavelength  and 
a  Density  Modulation  h  =  0.  4 
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